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Abstract 

The uniformly frustrated layered XY model is analyzed in its Villain form. A decouple 
pancake vortex liquid phase is identified. It is bounded by both first-order and second- 
order decoupling lines in the magnetic field versus temperature plane. These transitions, 
respectively, can account for the flux-lattice melting and for the flux-lattice depinning 
observed in the mixed phase of clean high-temperature superconductors. 

PACS Indices: 74.60.-w, 74.25. Dw, 74.25.Ha, 74.60.Ge 



The high-T c superconductor Bi 2 Sr 2 CaCu2 08 is perhaps the best known example of 
an extremely type-II layered superconductor. 1 In the absence of bulk pinning, the vortex 
lattice that exists in the mixed phase of this material melts across a first-order line in the 
temperature-magnetic field plane for magnetic fields applied perpendicular to the layers. 2 
The first-order line, H± = H m (T), begins at the zero-field critical point, T c , but it ends 
strangely in the middle of the phase diagram. The depinning line T = Td p (H±), which 
marks the point at which the flux lattice depins itself through thermal excitations, appears 
to be unrelated to this melting line. 3 ' 4 Considerable theoretical work has accompanied 
such observations. For example, flux-lattice melting is observed in Monte Carlo (MC) 
simulations of the frustrated XY model, 5 as well as in simulations of Ginzburg-Landau 
theory. 6 The Josephson coupling between layers is also predicted to diminish substantially 
at perpendicular fields, B±, that exceed the dimensional (2D-3D) cross-over scale, 7 B]_ ~ 
$o/Ao- Here, A denotes the Josephson penetration length. Last, it has been claimed 
that the layers decouple completely through a first-order transition at perpendicular field 
components H±_ many times larger than B*[. 8 A theoretical explanation of the critical 
endpoint of the first-order melting line in clean high-T c superconductors is lacking, however. 

In this paper, we obtain a schematic phase diagram for the uniformly frustrated 
XY model composed of a finite number of weakly coupled layers, which can describe the 
mixed phase of extremely type-II layered superconductors. 5 ' 9-12 In contrast to previous 
work, 7 ' 8 ' 12 the duality analysis 13-16 that follows includes both gaussian and topological 
excitations 17-19 within the intra-layer vortex lattices. We find a second-order melting line 
for weakly coupled two-dimensional (2D) vortex lattices in the limit of high perpendicular 
field. 17 ' 20 It is then argued that this line ends at an intersection with a first-order decou- 
pling line 8 in the middle of the T-H± plane (see Fig. 1). The theory potentially accounts 
for the same phenomenon that is observed in clean high-T c superconductors. 2 ' 3 We also 
demonstrate how the decoupling mechanism explicit in the duality analysis is intimately 
related to the entanglement of flux lines perpendicular to the layers. This is consistent 
with the supersolid scenario, 19 with recent numerical studies of the anisotropic XY model 
with frustration, 21 and with very recent experimental determinations of the phase diagram 
in clean high-temperature superconductors. 22 

Consider the application of a magnetic field to a layered superconductor in the extreme 



type-II limit, 1 A^ — > oo. Magnetic screening effects are then negligible and the theoretical 
description of the interior of the mixed phase reduces to a layered XY model in the 
presence of a uniform frustration. 5 ' 9-12 In particular, the thermodynamics is determined 
by the kinetic energy Exy = — J2 r n ^ cos [^a'^ , — Ad as a function of the superconducting 
phase 4>(r). Here, J x = J» = J y and J z = J± are the local phase rigidities, and A M = 
(0, b±x, — 6||x) is the vector potential. The magnetic induction parallel and perpendicular to 
the layers is related to the frustration, b, through the respective identities B\\ = (3>o/2Trd)b\\ 
and B± = (<& /27ra)b±. Here a denotes the square lattice constant, which is of order the 
zero-temperature coherence length, while d represents the spacing in between consecutive 
layers. Consider now the corresponding partition function, Z\p\ = J T)(pe~ ExY ' kBT e t 2^ p ^ 7 
in the (dual) Villain form, 13-16 which reads 
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Above, n^(r) is an integer link- field on the layered lattice structure of points r = (r, /), with 
fj, — x,y,z and n = (n x ,n y ). Also, we set (3\^j_ = J|| 5 j_/A;_bT. Now decompose the parallel 
field ft into transverse and longitudinal parts n{f, I) = n'(r, 1) — n-(f, I) + n-(r, I — 1), 
where the transverse and longitudinal fields, n' and n_, respectively satisfy the constraints 
V • n' = and V • n_ = n z (for p = 0), with V = (A x , A y ). 12,23 Take next the potential 
representation n_ = — V$ for the longitudinal inter-layer field. This yields the expression 
$(r, /) = ^2^/ G ( 2 )(r — r')n z (r', I) for the potential, where G^> = —V -2 is the Greens 
function for the square lattice. We then obtain the form Z[0] = Zqg ' H;Zdg[0] for the 
partition function (1), where 
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is an inter-layer Coulomb gas (CG) factor. Here, qi(r) = n z (r, / — 1) — n z (r, I) is the fluxon 
charge that collects onto layer / and yo = exp(— 1/2(3±) is the fugacity that is raised to 
the power N[n z ] = ^fW^^O P er configuration. 16 Note that a fluxon, n z (r,l) = ±1, 
represents a vortex ring that lies in between layers / and / + 1 at the point r. 24 The 



remaining factors 
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above represent modified 2D XY models with uniform frustration corresponding to each 
layer, with a modified in-plane vector potential A' = A + i/37 [n_(r, / — 1) — n_(r, /)] in the 
gauge V- A = 0, and with an extra weight factor C sw [qi] = exp{— J2^P\\)~ 1 [n-(f, / — 1) — 
fi-(r, I)] 2 }- In physical terms, the CG factor (2) describes the Josephson coupling between 
layers, whereas the discrete gaussian (DG) model factors (3) describe the thermodynamics 
of pancake vortices within each layer. The last factor in expression (2), however, represents 
the renormalization of the Josephson coupling due to misalignments of pancake vortices 
between layers. 7 ' 8 This important correction was omitted without proper justification in 
previous work. 12 ' 23-25 

Consider now the weak-coupling limit, yo — > 0, in which case the n z fluxons are 
dilute. 24 The modified 2D XY model (3) can then be analyzed in the continuum limit, in 
which case we obtain the identification 



^dg[«]/^dg[0] = (exp\i^2qi(r)(j) vx (rj) 
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i.e., 



with the vortex component of the phase correlations within an isolated layer / 
V 2 (^>vx = 0. It is instructive to consider a single neutral pair of unit n z charges that lie in 
between layers /' and I' + 1, separated by f. The renormalization to the Josephson coupling 
that is encoded in the last factor of expression (2) then becomes the gauge-invariant product 
nz(Zp G [<ft]/ZDG[0]) = C[, (r)C' l T +1 (r) of the corresponding phase autocorrelation functions, 
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within isolated layers. In the asymptotic limit, r — ► oo, this function has a magnitude 
of the form |C/(f)| = fi f o( r o/ r o) ??sw ( r o/l^*l)' ?vx f° r \A ^ £vx, and a magnitude of the form 
\C[if)\ = godrl/r'v) 11 ^ exp(-|f |/Cvx) for \r\ > ^ vx . 13 " 15 Here, ?y sw = (2tv(3\\)- 1 and rj vx are, 
respectively, the spin- wave and the vortex components of the correlation exponent inside 



layer /, while £ vx is the corresponding phase correlation length. Also, the length ro = 
a vx /2 3 / 2 e 7B is set by the inter- vortex scale, a vx = (&o/B±) 1 / 2 , and by Euler's constant, 14 
7s, while 7"o = a/2 3 ' 2 e lE . We therefore obtain the effective layered CG ensemble 

Z CG =J2y N[nz]e ^\-lYl E Qi(fi)[v2DHro/\?i -r 2 |)-^ s [ t i g (ri,f 2 )] g/ (r2) 
{n z } I / n/r 2 
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that has been coarse grained up to the natural ultraviolet scale a vx . In particular, the sums 
above are restricted to a square sublattice with lattice constant a vx . This requires the 
introduction of an effective coarse-grained fugacity y = go(a vx /a) 2 yo- At relatively small 
separation \fi — r 2 | <C £ vx , the fluxons experience a pure Coulomb interaction (V^ = 0) 
set by the 2D correlation exponent r]2D = %w + fe- At large separations \f\ — f^l ^ 
£ vx , on the other hand, the fluxons experience a pure [j\2d = 0) confining interaction 
Ktring(^i' ^ 2 ) = l^ 1 ~~ ^2|/Cvx between those points f\ and f-2 in layer / that are connected 
by a string [see Eq. (3) and ref. 26]. 

We shall now see that the effective layered CG ensemble (6) can be employed to 
determine the macroscopic nature of the Josephson effect in the weak-coupling limit, y — ► 0. 
Let us first identify the macroscopic intra-layer phase rigidity, 27 Jn = ksT '/2irr]2D- The 
above Coulomb gas ensemble indicates that fluxons (n z = ±1) are in a plasma state at 
temperatures below the naive decoupling temperature /c^T* = AirJn when quasi long-range 
intra-layer phase correlations are present, 23 ' 24 £ vx = 00 and V^ = 0, while that they 
form a dilute gas of bound neutral pairs of size £ vx at high temperatures when short- 
range intra-layer phase correlations exist, £ vx < 00. We now quote the expression for the 
perpendicular phase rigidity of the XY model (2) in terms of fluxons: 15 

/^^"'([E^'Of^BT, (7) 

r,l 

where M denotes the number of links between layers. It is obtained directly from the 
duality transformation (1) and from the definition pj: = ;^p-(G con d/A/') for this quantity. 
Under periodic boundary conditions, the low-temperature (plasma) phase thus sustains a 
macroscopic Josephson effect (pj- 7^ 0) due the presence of free fluxons, whereas the high- 
temperature (dielectric) phase does not due to the absence of free fluxons (J2 n z = 0)- 



Since the ordering temperature of a single layer is typically much less than T* (see be- 
low) , we conclude that the only thermodynamic phases that are possible at weak coupling 
are a coupled superconductor at low temperatures and a decoupled "normal" state at 
high temperatures. This indicates that neither the Friedel scenario 24 ' 28 (decoupled super- 
conducting layers) nor the line-liquid state 9 ' 12 ' 23 (coupled normal layers) are likely to be 
thermodynamic states in the absence of disorder. 

It is also important to determine the size of the local Josephson coupling. Consider 
the macroscopically decoupled phase at weak coupling, where short-range intra-layer phase 
correlations are present: £ vx < oo. Inter-layer fluxon (n z ) pairs are then bound by a 
confining string. Comparison of the CG ensemble (6) with the layered XY model (1) 
yields Koshelev's formula 10 

( e **«,« + i) =* yo J d 2 rC^C; +1 (r>-* 6 ii7a 2 (8) 

for the local Josephson coupling (see refs. 15 and 16). Here, 0^ + i(f) = <f)(r, I + 1) — 
4>{r,l) — A z (r) is the gauge-invariant phase difference between consecutive layers, while 
Ciif) is the phase autocorrelator for layer / in isolation [i.e., replace VX — ► in Eq. (5)]. 
Eq. (8) implies that d 2 \n.Z[O]/d0±dBn = at B» = 0. We then have a null line tension for 
Josephson vortices, since s\\ = = e|||jj.=o- I n conclusion, we recover the previous result 
(7) that macroscopic Josephson coupling is absent in the weak-coupling limit if intra-layer 
phase correlations are short range. Next, observe that scaling considerations yield the form 
J d 2 rCiCl +1 = /o£vx f° r the integral on the right-hand side of Eq. (8), where /o is of order 
unity. 10 Substitution into Eq. (8) then yields Koshelev's formula 

(cos0 M+1 )^/ o (£vx/a) 2 y o (9) 

for the local Josephson coupling in such case (see ref. 16). [Note that scaling implies the 
functional form £ vx = a vx e(/3||) for the 2D correlation length .] 

We can now analyze the XY model (1) composed of a finite number of weakly-coupled 
layers with uniform frustration, B±. Consider first the weak-coupling limit, (cos^z+i) — ► 
0, which by Eq. (9) is reached at infinitely high perpendicular fields. It is well known 
that an isolated lattice of 2D vortices (3) melts at a temperature ksTrn — ^||/20, above 
which quasi long-range positional correlation in the vortex lattice is lost. 20 The transition 
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is driven by the unbinding of dislocation pairs and it is expected theoretically to be second- 
order. 17 We shall now make the plausible assumption that the nature of phase coherence 
in the 2D vortex lattice is locked to the nature of the positional correlations, such that 
£ vx diverges exponentially as the temperature approaches T m in the disordered phase. 
By the previous analysis, we then conclude that the layers show a macrsoscopic Josephson 
effect at low temperature T < T m signalled by a positive phase rigidity between layers 
(7), while that they are decoupled (pjr = 0) at high temperature T > T m ' . The former 
low-temperature phase is best described by 2D vortex lattices that display a Josephson 
effect. The decoupled high-temperature phase, on the other hand, corresponds to a liquid 
of intra-layer "pancake" vortices. ' 8 

Consider next the weak-coupling regime, (cos0^+i) <C 1, at high perpendicular fields 
B±_ ^> B\_ [see Eq. (9)]. Eq. (9) indicates that the selective high-temperature expansion 
breaks down ((cos0^+i) > 1) in the decoupled phase at a temperature T x set roughly 
by the identification of length scales Aq ~ £vx(^x)- We now observe that the layered XY 
model (1) without frustration can also be described by the Coulomb gas ensemble (6), 
but with the natural ultraviolet length scale replaced globally by a vx — ► a. By analogy 
with what is presently understood for the layered XY model without frustration, 29 we 
conclude that a second-order transition should take place in the weak-coupling regime at 
a temperature T m that lies inside of the dimensional crossover window T m < T < T x . 

And what happens as the local Josephson coupling (9) approaches unity, which can be 
achieved by lowering the perpendicular field? The CG ensemble (6) is screened in the low- 
temperature Josephson-coupled phase, T < T m , for small effective fugacity. This implies 
that no phase transition is possible as a function of field there. 13 ' 1 Nevertheless, Eq. (9) 
clearly indicates that the selective high-temperature expansion, yo — ► 0, breaks down at 
perpendicular fields below B±. In such case, a crossover into a flux-line lattice regime must 
therefore take place. 7 At high temperatures T > T x , on the other hand, the CG ensemble 
(6) is confining for small fugacity, y C 1. In particular, the string interaction (6) binds 
together dilute fluxon-antifluxon pairs into stable dipoles of dimension £ vx . In the limit of 
dense fluxons, y — ► 1, these dipoles disassociate, however. This is due to the ineffectiveness 
of the string when the distance, r s , between neighboring dipoles is small in comparison to 
the length of the string, £ vx . The system must therefore experience a (inverted) first-order 



phase transition into a screened CG above a critical coupling (cos0^+i).d, since there is 
no diverging length scale. Monte Carlo simulation of the layered XY model with uniform 
frustration indicates that (cos^+i)^ is constant and of order unity. 11 By Eq. (9), we 
therefore expect a first-order decoupling transition at a perpendicular field 

H D = (/o$o/a 2 )(^/(cos</> M+1 > D ) (10) 

of order (3\\ B*^ for high temperatures, T > T x . This results from the replacement yo — > |/9j_ 
(see ref. 16). The phenomenology 30 Jj_ = Ejq(T c — T)/T c for the Josephson energy in the 
vicinity of the zero- field transition at T c yields the dependence 2 ' 7 ' 8 Hd(T) = 7^" H c2 (T) 
for the "cosine" XY model, where H C 2(T) ~ (<&o/a 2 )(T c — T)/T is the mean-field per- 
pendicular upper-critical field, and where 72 ~ ((e l< ^> !+1 ) d I fo) 1 ^ 2 • (ksTc/Ejo) 1 ^ 2 is an 
effective anisotropy parameter. Similar results for the first-order decoupling transition 
(10) were obtained previously using the elastic medium description of vortex matter in 
layered superconductors. 7 ' 8 The above discussion is summarized by the schematic phase 
diagram in Fig. 1. 

We shall now give a physical interpretation of the results just obtained from the 
duality analysis. Consider two parallel vortices along the z-axis that exchange positions in 
between layers / and / + 1. A moment's thought determines that the exchange creates a 
(fluxon) vortex loop that lies in between those layers. Fluxon charge 24 (n z ) can therefore 
entangle vortex lines aligned perpendicular to the layers. We conclude that entanglement is 
what actually drives the decoupling transitions shown in Fig. 1. This picture is consistent 
with (a) the classification of the coupled 2D vortex-lattice phase as a type of super-solid 
matter, 19 with (b) recent Monte Carlo simulations of the anisotropic XY model with 
frustration that conclude that entanglement is what drives the first-order melting transition 
of the vortex lattice, 21 and with (c) recent experimental work that also finds entanglement 
to be what controls the location of the critical endpoint of the vortex-lattice melting 
transition in clean high-temperature superconductors. 22 This web of facts strongly supports 
the results obtained here. 

We shall close by comparing the present theory for the uniformly frustrated layered 
XY model with known experimental results for the mixed phase of E^S^CaC^Os, which 
is extremely type-II and layered. This system shows first-order melting in the absence of 
bulk pinning. In agreement with Fig. 1, the melting line ends in the middle of the phase 



diagram. 2 This suggest identifying it with the decoupling transition at H± = Hp. Also, the 
depinning line, T = T^p, is nearly vertical in such materials for perpendicular fields above 
S^. 3 ' 4 This suggests identifying Td p with T m in Fig. 1. Finally, although Monte Carlo 
simulations of the frustrated XY model with anisotropy do observe a unique first-order 
vortex-lattice melting transition, no indication of the critical endpoint predicted here has 
been reported. 11 Eq. (10) implies, however, that the critical anisotropy parameter, 7' = 
(Jy/JjJ 1 / 2 , is high: e.g., Y c ~ 33 for perpendicular fields B± = $ /56a 2 at T m . The 
condition that the Josephson penetration length, Aq = 7^0, be smaller than L/2tt, where 
L is the linear dimension of each layer, then indicates that extensive MC simulations 21 are 
necessary in order to see the critical endpoint predicted here. 15 
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Figure Captions 

Fig. 1. Shown is a schematic phase diagram for the uniformly frustrated XY model (1) made 
up of a finite number of weakly-coupled layers at B\\ = 0. A vestige of the vertical 
second-order line may extend down to lower fields in the form of a cross-over (compare 
with ref. 3). The Josephson temperature Tj = J±/kB is assumed to be smaller than 
the scale of the figure. Notice that the inequality Tj < Tm required by the phase 
diagram indicates a minimum anisotropy parameter 7' = (J\\/J±) 1 ^ 2 equal to about 
four, since k B T^ D) £ Jii/20. 
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